Restoring site percolation on a damaged square lattice 
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We study how to restore site percolation on a damaged square lattice with nearest neighbor (N^) 
interactions. Two strategies are suggested for a density x of destroyed sites by a random attack 
at Pc- In the first one, a density y of new sites are created with longer range interactions, either 
next nearest neighbor (N^) or next next nearest neighbor (N'*). In the second one, new longer range 
interactions or N* are created for a fraction v of the remaining (pc — x) sites in addition to their 
interactions. In both cases, the values of y and v are tuned in order to restore site percolation 
which then occurs at new percolation thresholds, respectively tts, 7r4, 1123 and 7124. Using Monte 
Carlo simulations the values of the pairs {y, 713}, {y, 714} and {d, ■7r23}, {v, 7i"24} are calculated for the 
whole range < 2; < pc(N^). Our schemes are applicable to all regular lattices. 

PACS numbers: 07.05.Tp, 64.60.Ak 
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For several decades the calculation of percolation 
thresholds has been an ongoing technical challenge (see 
for the most recent ones). Up to now, analytical so- 
lutions have been limited to few a very two-dimensional 
lattices. However computer Monte Carlo simulations 
have allow to estimate percolation thresholds for an in- 
creasing large number of systems H, ^ IE 111 - High tem- 
peratures series expansions have been also extensively 
used 'T]. The most studied case is the hypercube due 
to its simple computer representation with available nu- 
merical estimates of percolation threshold up to thirteen 
dimensions 'sj, ^ 0- Some universal formulas have been 
also suggested 13, E|- The percolation phenomenon 0,01 
may also occur on a large variety of complex networks 
where sites are not distributed regularly [ij [13, ■ The 
site coordination number z varies for each site with a dis- 
tribution which depends mainly on the network growth 
rules. For example, when subsequent sites are attached 
randomly to alre ady existing ones the site degree becomes 
exponential [13 . H^ j . When the attachment is preferen- 
tial 0, ^3 the distribution follows a power law, and for 
so-called classical random graphs it is given by a Pois- 
son distribution |l2l Il6| . Such growing networks may 
reflect some properties of real networks — like WWW or 
Internet infrastructure — and have been the subject of 
studies of their resistance on intentional or random at- 
tacks At odd and for the first time we investigate 
in this paper the vulnerability at pc of a regular lattice 
with nearest neighbor (N^) interactions under a random 
attack. In particular we study how to restore site perco- 
lation on a square lattice once a fraction x of the initial 
Pc sites have destroyed. Two strategies are suggested. In 
the first one, new sites are created with longer range in- 
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FIG. 1: Various sites neighborhoods on the square lattice: (a) 
N^, (b) N^, (c) N*, and examples of their combinations: (d) 
N2-hN^ (e) N^+N'', (f) N^+N-^+N*. 



teractions, either next nearest neighbor (N'^) or next next 
nearest neighbor (N^). It is worth to note that when we 
consider sites with N'^ or N'' interactions, only N^^ or 
interactions are active. It means in particular that for 
these sites there exist no interactions. 

In the second strategy, no additional sites are created 
but instead new longer range interactions are created for 
a fraction v of the remaining not damaged after attack 
{pc — x) sites. We consider N"^ or interactions in addi- 
tion to their former interactions. Therefore v{pc — x) 
sites have plus either N'^ or interactions while 
{\ — v){pc — x) sites have only their initial interactions. 
Fig. nshows sites with N^, N^, N^, (N^+N^), (N^-l-N'') 
and (N^-l-N'^-l-N^) neighborhoods on the lattice. 

In both cases, the values of y and v are tuned in or- 
der to restore site percolation which occurs then at new 
percolation thresholds, respectively 7r3, 7r4, 7723 and 7r24. 
Using Monte Carlo simulations the values of the pairs 
{j/i'^'a}, {yi''^4} a-iid {i',7r23}, {w,7r24} are calculated for 
the whole range < a; < pc(N^). Our both schemes are 
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FIG. 2: (a) The percolation threshold ttc dependence on the 
neighborhoods mixing parameter q. In part (b) the results for 
TTa = 7rc(g,N^,N^) and n4 = 7rc(g, N^, N*) are repeated. The 
lines are guides for eyes. 



applicable to all regular lattices. 

However, as it will be explained below, it is much eas- 
ier, for computational reasons, to evaluate the percola- 
tion threshold on a square lattice with a given fraction 
q of occupied sites with one kind of neighborhood while 
remaining fraction (1 — 5) of occupied sites have another 
one. Then, having the percolation threshold tTc and the 
mixing neighborhood parameter q we can easily extract 
thevalues of {x,y} or {x,v}. We also estimate percola- 
tion thresholds tTc for the different values of the mixing 
parameter q for various pairs of mixed neighborhoods, re- 
spectively with N^, with N^, with (N^-^N^), 
with (N^-HN*), (N2-hN3) with (N^-hN"*), (N^-^N^) with 
(N2+N3-hN4) and (N^+N^) with (N^+N^+N^). These 
thresholds are presented in Fig. |21 

To perform our calculations of percolation thresholds 
we are using the Hoslicn-Kopclman algorithm (HKA) 
[TtI I from existing computational techniques 18, 19, 20]. 
With HKA each occupied site gets a label. The sites 
in the same cluster have the same labels and different 
labels are assigned to different clusters thus allowing to 
recognize which sites belong to which clusters. The HKA 
becomes particularly efhcient for the bus-bar percolation 



FIG. 3: Are full-filled black and red sites in one cluster? The 
red site has while the black one N^ Our answer is posi- 
tive, while the HKA will give such an answer only with the 
probability of 50%. The answer is sites labeling order depen- 
dent. 



problem, i.e. when we check if the site in distance i from 
the first fully occupied line is still connected to that line 
through the sites at the distances smaller than £. Then 
a computer stores only single line of sites and may go 
through the lattice only once . However, when all pos- 
sible links between sites are desired we have to keep the 
whole lattice during the simulation 0, |^ . The situation 
becomes even more complicated when we mix neighbor- 
hoods on the lattice, i.e. when different sites have differ- 
ent neighborhoods, as presented in Fig. |21 There, black 
site is not of dark (red) one, while dark (red) one is 
of black one. Note, that both sites are in the same clus- 
ter. However, the HKA will classify them as the mem- 
bers of one common cluster when it goes in type-writer 
order , and will leave them separated when it goes in 
reverse- type- writer order Thus, for the lattices with 
mixed neighborhoods we consider each site as having ef- 
fectively only slashed (red) 'half of their neighborhood 
presented in Fig. ^ due to computational reasons only. 

To start our simulations we denote 7rc(g,Ni,N2) the 
percolation threshold for a lattice on which the fraction 
q of occupied sites has Ni neighborhood while the re- 
maining fraction (1 — g) of occupied sites has N2. The 
percolation thresholds values tTc are then evaluated from 
the crossing point of two curves showing dependences of 
the percolation probability P vs. the sites occupation 
probability p for lattices of linear sizes L = 100 and 500 
(see Fig. 0)). The results are averaged over N^^n = 10'^ 
and 10^ for L = 500 and 100, respectively. They are pre- 
sented in Fig. 12 As a matter of fact, each site has only 
'half of their neighborhood as we have explained earlier. 



At this stage we are in position to determine which 
site density y with long-range interactions (for instance 
N3 qj. ]\[4^ necessary to reconstruct a percolation phe- 
nomenon which occurs now at = 7rc((7, N^, N"^) or 
774 = 7rc((7, N^, N"'), when the fraction x of the occupied 
sites with the neighborhoods has been removed due 
to the random attack. The results for selected values 
of parameter q are presented in Fig. [2Jb). They may 
be easily transformed to answer that question for both 
and additional interactions. To do that, we esti- 
mate the percolation threshold tTc for given pair of neigh- 
borhoods (here with and with N*) and given 
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7iJq,N^,N^+N'), qfrom 0.1 to 0.9 every 0.1, L=100 (symbols), L=500 (lines) 
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FIG. 4: The percolation probability P dependence on the 
occupation probability p for different values of mixing pa- 
rameter q which changes from 10% to 90% every 10% from 
right to left. The N'^ and (N^+N^) neighborhoods are mixed. 
The symbols correspond to L = 100 while lines to i = 500. 
The crossing points predict the percolation thresholds 7123 — 
7^,(5,N^N2 + N3). 



TABLE I: The fraction y of sites with which must be add 
to reconstruct the percolation phenomenon — which occurs 
at TTa = TVc{q, N^, N"^) — when only [pc(N^) — x] of sites with 
is occupied. 



1-q 


X 


y 


7i-3 


0.0 


0.0000 


0.0000 


0.592 


0.1 


0.0574 


0.0594 


0.594 


0.3 


0.1783 


0.1773 


0.591 


0.4 


0.2404 


0.2344 


0.586 


0.5 


0.3025 


0.2895 


0.579 


0.6 


0.3640 


0.3420 


0.570 


0.8 


0.4810 


0.4440 


0.555 


0.95 


0.5645 


0.5254 


0.553 


0.99 


0.5863 


0.5603 


0.566 


1.0 


0.5920 


0.5920 


0.592 



ratio q/{l — q) of their homogeneous mixture on the lat- 
tice. For estimated value of tTc we find immediately the 
fraction y = (1 — q)7:c of sites with long-range interac- 
tions. The remaining fraction of occupied sites, which is 
qiTc = Pc(N^) — 2;, has kind of neighborhood. The 
results are collected in Tabs. and ^ and presented in 
Fig. Eta). 

The other strategy which may be applied for the re- 
construction process is to enrich the fraction v of sites of 
those which survived an attack by additional long-range 
bonds to other existing and not damaged sites. Again, 
numerically we deal with the square lattice on which a 
given fraction q of occupied sites has neighborhoods, 
while remaining (1-g) has N^-l-N^ or N^-l-N"*. The latter 
number is also the fraction v of sites with long interac- 
tion among all occupied sites. If the percolation for such 



TABLE IL The fraction y of sites with N** which must be add 
to reconstruct the percolation phenomenon — which occurs 
at 714 = 7rc(q,N^,N^) — when only [pc(N^) — x] of sites with 
is occupied. 



1-g 


X 


y 


71-4 


0.0 


0.0000 


0.0000 


0.592 


0.05 


0.0401 


0.0291 


0.581 


0.2 


0.1512 


0.1102 


0.551 


0.3 


0.2189 


0.1599 


0.533 


0.4 


0.2812 


0.2072 


0.518 


0.5 


0.3395 


0.2525 


0.505 


0.7 


0.4456 


0.3416 


0.488 


0.9 


0.5424 


0.4464 


0.496 


0.99 


0.5866 


0.5376 


0.543 


1.0 


0.5920 


0.5920 


0.592 



TABLE IIL The fraction v of sites with which must be 
enriched with or N* bonds to reconstruct the percolation 
phenomenon — which occurs at respectively 7r23 or 7r24 — 
when only [pc(N^) —x] of sites with N'^ was saved after attack. 



1-g 


2:23 


7I"23 


X24 


7r24 


V 


0.0 


0.000 


0.592 


0.000 


0.592 


0.0 


0.1 


0.017 


0.575 


0.033 


0.559 


0.1 


0.2 


0.035 


0.557 


0.065 


0.527 


0.2 


0.3 


0.054 


0.538 


0.097 


0.495 


0.3 


0.4 


0.073 


0.519 


0.126 


0.466 


0.4 


0.5 


0.092 


0.500 


0.153 


0.439 


0.5 


0.6 


0.110 


0.482 


0.178 


0.414 


0.6 


0.7 


0.130 


0.462 


0.200 


0.392 


0.7 


0.8 


0.149 


0.443 


0.220 


0.372 


0.8 


0.9 


0.167 


0.425 


0.239 


0.353 


0.9 


1.0 


0.185 


0.407 


0.255 


0.337 


1.0 



mixture occurs at ttc then all of these sites have in 
their neighborhoods. The fraction x which is missed to 
the Pc(N^) is then simply pc(N^) — tTc- The results are 
collected in Tab. Illll and presented in Fig. IHb). 

For 'pure' cases, when each site has the same neigh- 
borhoods (i.e. when q = or q = 1) the results agree 
very well with thresholds Pc recently presented in Ref. 
I23I Similar checks for bi-colored lattices with and 
N'^ neighborhoods were performed in Ref. |3 For sites 
reconstruction process, when a: = or x = Pc(N^), the 
difference between originally used sites and their 
and/or substitutes vanishes: 

• In the limit a; — > the damages are very rare and 
thus reconstruction process is not too costly inde- 
pendently what kind of 'bricks' (sites with or 
N*) is used. 

• On the other hand, when x Pc(N^) the system 
must be completely reconstructed starting from 
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FIG. 5; The fraction of sites which must be (a) reoccupied 
(y) or (b) enriched with long-range interactions (v) to recover 
percolation phenomenon when the fraction x of sites with the 
was emptied. The lines are guides for eyes. 



'zero'. Since all of 'pure' N^, and N"* lattices has 
the same percolation threshold Pc(N^) = Pc(N^) — 
p^(N'*) = 0.5927460 ^ the cost of reaching 
the percolation threshold becomes again neighbor- 
hoods independent. 

• Note, that in both cases (N'^ and N"*) there is an 
optimal ratio q/{l — q) for which the percolation 
threshold tTc is the lowest. 

• For sites used for the reconstruction process 
^3 > Pc(N^) when 0.8 < g < 1, i.e. when dam- 
ages are relatively small (x < 0.2). In that case 
addition sites with diagonal bonds do not much 
help with restoring the percolation but increase to- 



tal sites density. 

For the second scenario, there is a critical density of 
sites Xc above which damages are not compensable by 
adding long range bonds even to all of remaining sites. 
For X > Xc the sites concentration falls below the thresh- 
old Pc for the square lattice with all sites having the same 
mixed neighborhood, for instance N^-l-N'^ or N^-l-N"* l23j . 
The critical values Xc(N) = 7rc(0,N^,N^-hN)-pc(N^) are 
0.185, 0.255 and 0.304 for N=N3, and N^-hN^, respec- 
tively. 

For both reconstruction strategies, the percolation oc- 
curs more easily, it is achieved at a lower cost, if neigh- 
borhood with larger radius are employed as presented in 
Fig. El 

To conclude, in contrast to Ref. we have not been 
interested on how to destroy connectivity on a given lat- 
tice but rather on how to effectively reconstruct it after a 
random destruction. Using Monte Carlo simulations we 
have reported square lattice site percolation thresholds 
TTc for given neighborhoods characterized by a mixing pa- 
rameter q and various pairs of mixed neighborhoods build 
from basics ones, i.e. N^, and N''. We have showed 
quantitatively that 'repairing costs' are lower when sites 
with longer radius of neighborhood are employed for site 
reconstruction process at intermediate range of damages. 
The strategy involving bond enrichment fails if damages 
are too large. The critical damages size Xc depends on 
the percolation threshold Pc(N^ -f N), where N is the kind 
of bonds used in the enriching process. Both our schemes 
are readily applicable to all regular lattices. The re- 
sults may found useful applications in the fields of econo- 
physics 0] and sociophysics 0, for instance to study 
reconstruction processes in social systems after an expo- 
sion to some sort of damages where interpersonal bonds 
have been broken. 
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